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Abstract 

Bose-Einstein correlations have been measured using samples of proton-proton colli- 
sions at 0.9 and 2.36 TeV center-of-mass energies, recorded by the CMS experiment at 
the CERN Large Hadron Collider. The signal is observed in the form of an enhance- 
ment of pairs of same-sign charged particles with small relative four-momentum. The 
size of the correlated particle emission region is seen to increase significantly with the 
particle multiplicity of the event. 



*See AppendixfAJfor the list of collaboration members 



In particle collisions, the space-time structure of the hadronization source can be studied using 
measurements of Bose-Einstein correlations (BEC) between pairs of identical bosons. Since the 
first observation of BEC fifty years ago in proton-antiproton interactions [IJ, a number of mea- 
surements have been made by several experiments using different initial states; a detailed list 
of the experimental results can be found in [21131 ■ Boson interferometry at the Large Hadron 
Collider provides a powerful tool to investigate the space-time structure of the particle emis- 
sion source on femtometric length scales at different center-of-mass energies and with different 
initial states, using the same detector. This letter reports the first measurement of BEC param- 
eters in pp collisions at 0.9 and 2.36 TeV w^ith the CMS detector. 

Constructive interference affects the joint probability for the emission of a pair of identical 
bosons with four-momenta pi and p2- Experimentally, the proximity in phase space between 
final-state particles is quantified by the Lorentz-invariant quantity Q = -\/— (pi — ^2)^ = 
yM-^ — 4m^, where M is the invariant mass of the two particles, assumed to be pions with 
mass niji- The BEC effect is observed as an enhancement at low^ Q of the ratio of the Q distribu- 
tions for pairs of identical particles in the same event, and for pairs of particles in a reference 
sample that by construction is expected to include no BEC effect: 

R{Q) = (dN/dQ)/(dN,ef/dQ), (1) 

which is then fitted with the parameterization 

R(Q)=C[l + AQ(Qr)]-(l + ^Q). (2) 

In a static model of particle sources, Q ( Qr ) is the Fourier transform of the spatial distribution of 
the emission region of bosons with overlapping w^ave functions, characterized by an effective 
size r. It is often parameterized as an exponential function, Q(Qr) = e^^'', or with a Gaussian 
form, Q(Qr) = e^^'^''' (IHJ and references therein). The parameter A reflects the BEC strength 
for incoherent boson emission from independent sources, S accounts for long-range momentum 
correlations, and C is a normalization factor. 

The data used for the present analysis were collected by the CMS experiment in December 
2009 from proton-proton collisions at center-of-mass energies of 0.9 and 2.36 TeV. A detailed 
description of the CMS detector can be found in |5|. The central feature of the CMS appa- 
ratus is a superconducting solenoid of 6 m internal diameter, providing a uniform magnetic 
field of 3.8 T. The inner tracking system is the most relevant detector for the present analy- 
sis. It is composed of a pixel detector with three barrel layers at radii between 4.4 and 10.2 cm 
and a silicon strip tracker with 10 barrel detection layers extending outwards to a radius of 
1.1 m. Each system is completed by two endcaps, extending the acceptance up to a pseudo- 
rapidity I?/ 1 = 2.5. The transverse-momentum (pj) resolution, for IGeV charged particles, is 
between 0.7% at tj = and 2% at |//| = 2.5. The events were selected by requiring activity in 
both beam scintillator counters 16J. A minimum-bias Monte Carlo (MC) sample was generated 
using PYTHIA (with D6T tune) IT] follow^ed by full detector simulation based on the Geant4 
program |8|. Additional PYTHIA MC samples w^ere generated to simulate BEC effects with 
both Gaussian and exponential forms of Q(Qr). 

Charged particles are required to have px > 200 MeV, which is sufficient for particles emitted 
from the interaction region to cross all three barrel layers of the pixel detector and ensure good 
two-track separation. Their pseudorapidity is required to satisfy | /^track I < 2.4. To ensure high 
purity of the primary track selection, the trajectories are required to be reconstructed in fits 
with more than five degrees of freedom (dof) and X'^/^doi < 5.0. The transverse impact pa- 
rameter with respect to the collision point is required to satisfy \dxy\ < 0.15 cm. The innermost 



measured point of the track must be less than 20 cm from the beam axis , in order to reduce 
electrons and positrons from photon conversions in the detector material and secondary par- 
ticles from the decay of long-lived hadrons {K^,A, etc.). In a total of 270472 (13548) events 
selected at 0.9 (2.36) TeV center-of-mass energy, 2 903 754 (188 140) tracks are accepted by these 
selection criteria. 

All pairs of same-charge particles with Q between 0.02 and 2 GeV are used for the measure- 
ment. The lower limit is chosen to avoid cases of tracks that are duplicated or not w^ell sepa- 
rated, while the upper limit extends far enough beyond the signal region to verify a good match 
between signal and reference samples. A study with simulated data shows that the ratio of the 
tracking efficiencies of particle pairs in the signal and in the reference samples is independent 
of Q in the measurement region. 

Coulomb interactions between charged particles modify their relative momentum distribution. 
This effect, which differs for pairs with same charge (repulsion) and opposite charge (attrac- 
tion), is corrected for by using Gamow factors [9J. As a cross-check, the enhancement in the 
production of opposite-charge particle pairs with small values of Q is measured in the data and 
is found to be reproduced by the Gamow factors to within ±15%. 

Different methods are designed to pair uncorrelated charged particles and to define reference 
samples used to extract the distribution in the denominator of Eq. ([l]|. Opposite-charge pairs: 
this data set is a natural choice but contains resonances (//, p, ...) which are not present in 
the same-charge combinations. Opposite-hemisphere pairs: tracks are paired after inverting in 
space the three-momentum of one of the two particles: (£, p) — ?► (£, — p) ; this procedure is ap- 
plied to pairs with same and opposite charges. Rotated particles: particle pairs are constructed 
after inverting the x and y components of the three-momentum of one of the two particles: 
{px,Py,Pz) —> {—px,—Py,Pz)- Pairs from mixed events: particles from different events are 
combined w^ith the following methods: i) events are mixed at random; ii) events w^ith simi- 
lar charged particle multiplicity in the same rj regions are selected; iii) events with an invariant 
mass of all charged particles similar to that of the signal are used to form the pairs. 

As an example, the ratios R{Q) obtained with the opposite-hemisphere, same-charge reference 
samples are shown in Fig. Illboth for data and simulation without BEC. A significant excess at 
small values of Q is observed in the data. Additional details are given in [lOJ. 
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Figure 1: Ratios R{Q) obtained with the opposite-hemisphere, same-charge reference samples 
for data (dots) and MC with no BEC effect (crosses). 



In order to reduce the bias due to the construction of the reference samples, a double ratio TZ is 

defined: 

R / dN/dQ \ / / dNuc/dQ \ 



nQ) 



Rmc \dN,,i/dQj / \dNMC,rei/dQj' 



(3) 



Table 1: Results of fits to the double ratios TZ for several reference samples, using the parame- 
terization of Eq. (|2| with the exponential form, for 0.9 TeV data (top) and 2.36 TeV data (bottom). 
Errors are statistical only, and quoted as if independent. 



Results of fits to 0.9 TeV data 


Reference sample 


p value (%) 


C 


A 


r(fm) 


<5(10-3GeV-^) 


Opposite charge 


21.9 


0.988 ± 0.003 


0.56 ± 0.03 


1.46 ±0.06 


-4 ±2 


Opposite hem. same ch. 


7.3 


0.978 ± 0.003 


0.63 ± 0.03 


1.50 ±0.06 


11 ±2 


Opposite hem. opp. ch. 


11.9 


0.975 ± 0.003 


0.59 ± 0.03 


1.42 ±0.06 


13 ±2 


Rotated 


0.02 


0.929 ± 0.003 


0.68 ± 0.02 


1.29 ±0.04 


58 ±3 


Mixed evts. (random) 


1.9 


1.014 ±0.002 


0.62 ± 0.04 


1.85 ±0.09 


-20 ±2 


Mixed evts. (same mult.) 


12.2 


0.981 ± 0.002 


0.66 ± 0.03 


1.72 ±0.06 


11 ±2 


Mixed evts. (same mass) 


17.0 


0.976 ± 0.002 


0.60 ± 0.03 


1.59 ±0.06 


14 ±2 


Combined 


2.9 


0.984 ± 0.002 


0.63 ± 0.02 


1.59 ± 0.05 


8±2 




Results of fits to 2.36 TeV data 


Reference sample 


p value (%) 


C 


A 


r(fm) 


<5(10-3GeV-i) 


Opposite charge 


57 


1.004 ± 0.008 


0.53 ± 0.08 


1.65 ±0.23 


-16±6 


Opposite hem. same ch. 


42 


0.977 ±0.006 


0.68 ±0.11 


1.95 ±0.24 


15 ±5 


Opposite hem. opp. ch. 


46 


0.969 ± 0.005 


0.70 ±0.11 


2.02 ± 0.23 


24 ±5 


Rotated 


42 


0.933 ± 0.007 


0.61 ± 0.07 


1.49 ±0.15 


58 ±6 


Mixed evts. (random) 


23 


1.041 ± 0.005 


0.74 ±0.15 


2.78 ± 0.36 


-40 ±4 


Mixed evts. (same mult.) 


35 


0.974 ± 0.005 


0.63 ± 0.10 


2.01 ± 0.23 


20 ±5 


Mixed evts. (same mass) 


73 


0.964 ± 0.005 


0.73 ±0.11 


2.18 ± 0.23 


28 ±5 


Combined 


89 


0.981 ± 0.005 


0.66 ± 0.07 


1.99 ±0.18 


13 ±4 



where the subscripts "MC" and "MC,ref " refer to the corresponding distributions from the MC 
simulated data generated without BEC effects. 

The results of fits of TZ{Q) based on the parameterization of Eq. pi with Q(Qr) = e^Q'' are 
given in Table [ij both for 0.9 and 2.36 TeV data. In the case of the opposite-charge sample, it is 
found that the region with 0.6 < Q < 0.9 GeV, containing a sizeable contribution of pairs from 
p — > 7r+ n^ decays, is not well described by the MC |ITO| . This region is therefore excluded from 
the fits w^ith this reference sample and also w^ith the combined sample defined below. 

As a cross-check, the dE/dx [11] measurements of particles in the tracker are used to select a 
sample enriched in nn pairs, and another sample w^ith one of the particles not consistent with 
the pion hypothesis. Figure |2]presents the double ratios for these two samples at i/s = 0.9 TeV, 
showing that an enhancement at small Q values is observed only in the case of identified nn 
pairs. 

As none of the definitions of the reference samples is preferable a priori, an additional, "com- 
bined" double ratio 7?.'^°"^^ is formed, where the data and MC distributions are obtained by 
summing the Q distributions of the seven corresponding reference samples. 

The distributions of 7^'^°"^'^ for 0.9 and 2.36 TeV data are shown in Fig. pi and the values of the fit 
parameters are given in Table IT] A large correlation is found between the parameters A and r, as 
w^ell as betw^een 5 and C (correlation coefficients of 0.82 and —0.97 at 0.9 TeV, respectively). The 
data are described by Eq. (J2| with an exponential form for Q(Qr), as shown by the solid lines 
in Fig. Island confirmed by the fit probability (p value) in Table llj The fit with a Gaussian form, 
Q(Qr) = e-(Q'')', which yields A = 0.32 ± 0.01, r = 0.98 ± 0.03 fm, does not correctly describe 
the 7^(Q) distribution, as shown by the dashed lines in Fig.p^and by a p value of lO^'^^. Gaus- 
sian shape fits also proved to offer a poor description of the data in previous measurements 
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Figure 2: Double ratios 7^(Q) for the 0.9 TeV data, using the opposite-hemisphere, same-charge 
reference samples for combinations enriched, using a dE/dx measurement, in pion-pion pairs 
(dots) and in pion-non-pion pairs (open circles), respectively. 

Although the values of r obtained in the exponential fits cannot be compared directly with 
results obtained with a Gaussian function, it should be noted for comparison purposes that the 
first moment of the Q(Qr) distribution corresponds to 1/r for an exponential shape and to -^ 
for a Gaussian form. Alternative functions, as defined in IIT3l[T5l[T6| , also describe the data w^ell 
w^ith similar p values. In particular for the Levy parameterization, Q(Qr) = e^^^''' , the fitted 
values are A = 0.93 ± 0.11, r = 2.46 ± 0.38 fm, and a = 0.76 ± 0.06, with a p value of 12.8%. 
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Figure 3: Fits to the double ratios TZ'^°"^^{Q) with exponential (solid lines) and Gaussian 
(dashed lines) functions, for 0.9 TeV (top) and 2.36 TeV (bottom) data. The range 0.6 < Q < 
0.9 GeV is excluded from the fits. 



The leading source of systematic uncertainty on the measurements arises from the fact that 
none of the reference samples is expected to give a perfect description of the Q distribution 
in the absence of BEC, and that none of them can be preferred or discarded a priori. The cor- 
responding contribution to the systematic error is computed as the r.m.s. spread between the 
results obtained for the different samples, i.e., ±7% for A and ±12% for r. The systematic un- 
certainty related to the Coulomb corrections is computed by propagating the measured ±15% 
agreement margin, resulting in ±2.8% variation for A and ±0.8% for r. The presence of a pos- 
sible bias introduced by the track reconstruction and selection requirements was studied by 
comparing the results obtained at the generator and reconstruction levels in the MC simulation 
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Figure 4: Values of the A (top) and r (bottom) parameters as a function of the charged-particle 
multiplicity in the event for combined (dots) and opposite-hemisphere, same-charge (open cir- 
cles) reference samples, at 0.9 TeV. The errors shown are statistical only. The points are placed 
on the horizontal scale at the average of the multiplicity distribution in the corresponding bin. 



that incorporates BEC effects. The differences in the fitted parameter values for the different 
reference samples are smaller than the statistical errors and no systematic bias is observed for 
r. No correction is therefore applied and no additional systematic error is included. For the 
2.36 TeV data the same relative systematic uncertainties as for the 0.9 TeV results are used, in 
view of the reduced size of the sample and the larger statistical uncertainties of the fit results. 

The BEC parameters measured with the combined reference sample are 

A = 0.625 ± 0.021 (stat.) ± 0.046 (syst.) and r = 1.59 ± 0.05 (stat.) ± 0.19 (syst.) fm at 0.9 TeV; 

A = 0.663 ± 0.073 (stat.) ± 0.048 (syst.) and r = 1.99 ± 0.18 (stat.) ± 0.24 (syst.) fm at 2.36 TeV. 

The possible dependence of the BEC signal on various track and event observables has been 
studied. A significant dependence of r on the charged-particle multiplicity in the event is ob- 
served for all reference samples. Here, the only mixed-event reference sample used is the one 
constructed by combining charged particles from events in the same multiplicity range. The 
fit parameters for the combined reference sample are given in Table |2] and shown in Fig. |4] as 
a function of the track multiplicity for the 0.9 TeV data. As an example, the results for the 
opposite-hemisphere, same-charge reference sample are also shown in Fig. |4J The systematic 
errors on A and r in each multiplicity bin are taken as the r.m.s. spread of the results obtained 
with the various reference samples. Due to the limited sample size of the 2.36 TeV data only 
two multiplicity bins are considered, one for multiplicities smaller than 20 tracks, the other for 
multiplicities between 20 and 60 tracks. The values measured for the parameters with the com- 
bined reference samples are A = 0.65 ± 0.08 and A = 0.85 ± 0.17, and r = 1.19 ± 0.17 fm and 
r = 2.85 ± 0.38 fm for these two multiplicity bins, where the errors are statistical only. For com- 
parison, the values obtained for the same multiplicity bins at 0.9 TeV are A = 0.65 it 0.02 and 
A = 0.63 ± 0.05, and r = 1.25 ± 0.05 fm and r = 2.27 ± 0.12 fm, respectively. These measure- 
ments are consistent w^ithin errors. The dependence of r on multiplicity was already observed 
in previous measurements as discussed in detail in O. 

In summary, Bose-Einstein correlations have been measured for the first time at the LHC by the 
CMS experiment in pp collisions at 0.9 and 2.36 TeV center-of-mass energies. Several reference 
samples were used to extract the signal. For all of them an exponential shape fits the data 



Table 2: Results of the fits to the double ratio 'JZ'^°"^^ for the combined reference samples, using 
the parameterization of Eq. ^ with the exponential form, as a function of the charged-particle 
multiplicity in the event, for 0.9 TeV data. Errors are statistical only, except for A and r where 
statistical (first error) and systematic uncertainties (second error) are given. 



Mult. 


p val. 


C 


A 


r(fm) 


3 (10-^ 


range 


(%) 








GeV-i) 


2-9 


97 


0.90±0.01 


0.89±0.05±0.20 


1.00±0.07±0.05 


72±12 


10-14 


38 


0.97±0.01 


0.64±0.04±0.09 


1.28±0.08±0.09 


18± 5 


15-19 


27 


0.96±0.01 


0.60±0.04±0.10 


1.40±0.10±0.05 


28± 5 


20-29 


24 


0.99±0.01 


0.59±0.05±0.17 


1.98±0.14±0.45 


13± 3 


30-79 


28 


l.OOiO.Ol 


0.69±0.09±0.17 


2.76±0.25±0.44 


10± 3 



significantly better than a Gaussian shape. An increase of the effective size of the emission 
region with charged-particle multiplicity in the event has been observed. 
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